Generalized Einstein or Green-Kubo relations for active biomolecular transport 
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For driven Markovian dynamics on a network of (biomolecular) states, the generalized mobilities, 
i.e., the response of any current to changes in an external parameter, are expressed by an integral 
over an appropriate current-current correlation function and thus related to the generalized diffusion 
constants. As only input, a local detailed balance condition is required typically even valid for 
biomolecular systems operating deep in the non-equilibrium regime. 

PACS numbers: 87.16.-b, 05.40.-a 



Introduction.- Close to equilibrium, transport coeffi- 
cients like the mobility, conductivity or viscosity, quan- 
tifying the response of a system to an external field or 
perturbation, can be expressed by equilibrium correla- 
tion functions. The Stokes-Einstcin relation between the 
mobility and the diffusion constant of a spherical particle 
is arguably the oldest and best known example of such 
a Green-Kubo relation Both mobility and diffusion 
constant are still well-defined even for a non-equilibrium 
steady state (NESS) of an open or driven system in which 
stationary currents lead to permanent dissipation. In 
such a state, the Stokes-Einstein relation no longer holds 
true. The difference between diffusion constant and mo- 
bility, however, can be expressed by an integral over an 
experimentally measurable correlation function Q. 

In the present paper, we investigate the relation be- 
tween a mobility or transport coefficient and the cor- 
responding dispersion or fluctuations for any current in 
an arbitrary driven system with the special focus on 
biomolecular transport like the one mediated by molec- 
ular motors or ion channels and pumps. The essential 
characteristics of such transport is that even though the 
system is driven, typically by non-balanced chemical re- 
actions involving ATP, it takes place in a well-defined 
thermal environment. This fact imposes a constraint on 
the ratio between forward and backward rates for any 
mesoscopic transition that will allow us to express the 
difference between mobility and dispersion in a physically 
transparent way. On the technical level, we build on the 
recent derivation of a general fluctuation-dissipation the- 
orem for NESSs 0-0]. By directly working in the NESS, 
our approach is complementary to work that invokes the 
fluctuation theorem for deriving non-linear response co- 
efficients in higher order expansions around equilibrium 
0, [1]. Moreover, it goes beyond similar relations ob- 
tained for genuine diffusive spatial transport since 
we require no Euclidean metric and hence the notion of a 
locally co-moving frame is not available. Our results will 
therefore be applicable not only to any discrete model 
for a molecular motor or ion pump (see, e.g., [lll - flij and 
references therein) but also to driven (bio) chemical re- 
action networks and their response to changing chemical 
conditions 15, 3]. As a simple illustration will show, a 



misguided rewriting of our additive relationship between 
mobility and dispersion in terms of a multiplicative "ef- 
fective temperature" could easily lead even to negative 
values for the latter as found for various active biomolec- 
ular systems, see, e.g., [itI [isj. 

System.- We describe the system by a set of discrete 
states {n}. At time t, the system is in a state n{t) jump- 
ing at discrete times tj from state Uj" to state rij'^. A 
transition between state m and state n occurs with a rate 
Wmn- With each transition m n, we associate trans- 
port of a quantity = — (i"„j leading to a microscopic 
current 
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The transition rates between the states depend on a set 
of external parameters {/i^}. We make no particular as- 
sumptions on the parameter dependence of the individual 
transition rates but only require that the ratio between 
forward and backward rates obeys the typical "local de- 
tailed balance" (LDB) condition 
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which implies for the logarithmic derivatives, or "sensi- 



tivities", r^ = Tdh^ In Wmn, the crucial relation 
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Here, and throughout the paper, we set Boltzmann's con- 
stant /cb = 1. Examples for pairs of an external param- 
eter ha and a conjugate distance are (i) force / and 
spatial distance d, (ii) chemical potential and number 
of consumed (or, if negative, produced) molecules of 
type a (like ATP and ADP) and (iii) potential difference 
A(/) and transported electrical charge q. These choices 
are relevant to molecular motors (i-ii) and ion pumps (ii- 
iii), respectively. In all these cases, the LDB condition is 
usually assumed not only for small deviations from equi- 
librium but also for finite values of the fields {hp}. 

For constant external parameters {hp}, the system 
reaches a stationary state in which = {5n(t)m) de- 
notes the probability to find it in the particular state 



2 



m. Throughout the paper, the brackets (...) denote av- 
erages in this stationary state. If the system operates in 
a genuine NESS at least one pair of directed probabihty 
currents 



(4) 



is non-zero. Consequently, some of the currents have a 
non-zero mean 
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We will need a second type of current derived from a 
local variable 

''ait) = "^"(t)™^™ ^it'^ "^n = ^rnkr^ik/Pm 

m k 
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which could be called a "sensitivity- weighted" current. It 
generalizes the mean local velocity found in this context 
for Langevin systems Q to arbitrary networks. Positive 
contributions to Vait) arise from links for which the di- 
rected probability current and the sensitivity have the 
same sign. The dimension of i^a{t) justifies to call it a 
current. Moreover, its mean is equal to the ordinary cur- 
rent since (vait)) = Yjmn^rnnr",ri = 3a where we use 
([3]) from above. 

Generalized Green-Kubo relations.- The aim of gen- 
eralizing the Einstein or Green-Kubo relations to non- 
equilibrium processes requires that we express both the 
generalized diffusion constants, or dispersions, and the 
generalized mobilities by correlation functions involving 
currents. The dispersions given by 

lim 1 fdt' f dt"{(Ut')~]a){30{t")~Jp)) 
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characterize the integrated fluctuations around the mean 
currents. By isolating the diagonal in this double inte- 
gral and exploiting stationarity, we can rewrite the dis- 
persions as 

Dap = / dt{{ja{t) - j„)(jMO) - jp)) + D'°l. (8) 

The lower boundary 0+ at the integral indicates that no 
delta-like contributions at t = should be picked up since 
those are captured by the time-local contribution 
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The generalized mobilities Kafs = dh/jja quantify the 
dependence of the mean current on an external parame- 
ters. As our main result, we will prove below that they 



can also be expressed by an integral involving a correla- 
tion function of the currents just introduced and a local 
term in the form 
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Hence, the difference between the dispersion and mobility 
tensors can be expressed as 
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Jo 

(13) 

with the local contribution 
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indicates that each link is 



where the notation X]m<n 
counted only ones. 

In equilibrium, J^^(0),j/3 and Kmn a-U vanish identi- 
cally, and hence lap = 0. Our representation makes the 
"violation" of the Einstein or Green-Kubo in a NESS ap- 
parent and provides a physically transparent expression 
for the difference between dispersions and mobilities. 

Molecular motor - As an illustration of the general 
framework we consider any discrete state model of a 
molecular motor. A transition from state m to state n 
may either advance the motor a spatial distance dmn = 
—dnnn Or be associated with a chemical reaction of the 
type EaCn^a ^ Ea^mn^a, Or contaiu both. The 
index a = t,d,p labels the chemical species ATP, ADP 
and Pi, respectively, and and s^„(= r"^) are the 
corresponding stochiometric factors for the forward and 
backward reaction. For each transition and each species 
a "chemical distance" 
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denotes the number of consumed (or. if negative, pro- 
duced) molecules of type a. The chemical species are 
provided at externally controlled concentrations Ca- For 
any motor and no applied external force (/ = 0) there 
arc concentrations at which the motor is in equilib- 
rium with its thermal and chemical environment. As- 
suming ideal behaviour, the concentrations are linked to 
the chemical potentials by fia = p'^^ + T\xiCa/c°^- If, still 
at / = 0, the chemical potentials deviate from their equi- 
librium value, the transition rates are modified according 
to the usual mass action law kinetics. 
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where A/i^ = /i„ — Note that the dependence of 

these rates on the chemical potentials {ha = A/Zq) obeys 
^ which justifies a posteriori to denote the stochiometric 
coefficients by r^n- We make no particular assumptions 
on the force dependence of the individual transition rates 
but require that the ratio between forward and backward 
rates obeys, as usually assumed, the LDB condition 
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Hence, the sensitivities r„j„ = Tdf Inwmn obey the rela- 
tion (121). 

For a simple but still instructive specific example, we 
consider a "one-state" ratchet model where the forward 
rate (driven by ATP hydrolysis) and the backward rate 
(synthesizing ATP from ADP and Pi) are given by 
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and 
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w;"^exp[(A/ + AAi''-/6l_d)/r], (19) 



respectively. The load sharing factors 6+ and 6'_ with 
9+ + 9- = 1 guaranteeing the LDB condition ([3]) are 
related to the distance of the activation barrier in forward 
and backward direction, respectively [l2j . 

Since in this model all sites are physically equivalent 
but only spatially translated a distance d, there are no 
current correlations, so that only the local terms con- 
tribute. With j — d{'w^ — W-), the ordinary spatial 
mobility becomes /i = df j = d'^[9-i^w^ + 9^w^]/T, and 
the corresponding diffusion coefficient D = {l/2)d'^[w+ + 
W-]. The difference / = —d'^{w+ — w-){9^ — 9^)/2 van- 
ishes not only in equilibrium (?«+ = w_) but even in a 
NESS for a symmetric barrier (6*+ = 0_ = 1/2). 

Expressed in terms of an effective temperature, 
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where p = w+fw-, one sees that for < 6*+ < 1, T /T 
can acquire any value > 1/2. If we allow the somewhat 
more extreme structural choice of 0+ > 1 (thus assuming 
that both forward and backward steps are promoted with 
increasing force) then even negative values of the effective 
temperature become possible. Clearly, even this simple 
example demonstrates that the idea of phenomenologi- 
cally characterizing active processes by an elevated "ef- 
fective temperature" is not really consistent. It rather 
conceals the physically transparent additive relationship 
between mobility and dispersion by replacing it with a 
multiplicative factor. 

Rather than looking at the response of the motor to 
a changing applied force, one can ask for the response 
to a change in concentration of ATP or ADP, i.e., to a 
change in the chemical potential with hp = pp. For the 
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TABLE I: Generalized mobilities for the one-state motor. 



current, we can either choose the ordinary spatial current 
j{t) or the current of consumed a-moleculcs ja{t)- How 
the corresponding mean currents change with the chemi- 
cal potential of /3-molecules is expressed by the mobility 
tensor Kap shown in Table I which includes the "cross" 
mobilities between chemical and mechanical (here denote 
by an index /) distances and fields. We refrain from list- 
ing the dispersions, which are in this case symmetrical 
with Dai3 ~ Dpa, and the corresponding effective tem- 
peratures except for pointing out that the latter are 
asymmetric and depend on the choice of indices even for 
fixed rates. 

While the evaluation of mobilities and dispersions is 
straightforward also for any more complex specific model 
as will be illustrated elsewhere, a few universal state- 
ments seem to be possible beyond the obvious ones refer- 
ing to equilibrium. As one example consider the observa- 
tion made in [ll| for a particular two state motor model 
that at stalling conditions, j ~ at f = fg, the usual 
Einstein relation between mobility and diffusion constant 
holds true, even though idle chemical currents dissipate 
energy. Our expressions dH), ([13]) and ([14]) show that, 
in general, the validity of the Einstein relation requires 
not only that j = but moreover that any link carry- 
ing a non-zero probability current Kmn be not sensitive 
to the force /, i.e. for any K^n 7^ 0, = fnm = 
must hold. The latter condition will not necessarily be 
met at stalling since even pure chemical transitions with 
dmn = will, in general, be affected by changing the 
applied force. 

Proof of |77)) .- In the differential mobility 



Kq/3 = dhi^ia=^d'^^dhi,{PmWmn) (21) 
mn 

= Y.dZ,^{dh,Pm)Wmn + n'°l (22) 



the term dhfjPm = (i5jjYt)m) must be expressed by a 
correlation function. In [Gj we have determined the re- 
sponse of an observable at time t2 to a delta-like pertur- 
bation at time ti. Specialized to the present quantities 
and slightly adapting the notation, this relation reads 



t2)inl 



6hp{ti 



{/i^ }— const 



{5n{t,)mB{h)), (23) 



4 



where the conjugate variable B{ti) is given by 



where we have used (P) and ^ . If the correlation func- 



tion {Sn{t2)niYl,jHtl 



J ' ri r 



_) is averaged over the 



states n- before the jump at ti one gets 

{5n(t,)mY.^{tl-ty^+^-) (26) 
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on the external perturbation. While this is not an is- 
sue in any theoretical modelling, it will limit the direct 
application to those experimental systems for which this 
property of the rates is accessible. In the familiar linear 
response realm of the regular Green-Kubo relations, such 
explicite knowledge is not necessary. This observation 
might support the view that often the quantitative eval- 
uation of exact non-equilibrium relations requires more 
specific input than their equilibrium counter-parts do. 
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Putting together (|23ll27p and using (jS]), we can write 



T 



^hpi }— const 
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(28) 

Thus the response of the current at the later time, 
ja(^2) = J2ni^n{t2)ndniW„i, to a delta-like perturbation 
at the earlier time can be expressed as 



T 



{his} — const 



{jc.it2)Upiti)-Mti))- (29) 



Integrating over the time- difference t2—ti, we obtain our 
main result pT|) . 

Concluding perspective - We have expressed the gen- 
eralized mobilities by current correlation functions for 
any driven system described by a master equation with 
transition rates which obey a local detailed balance con- 
dition as it should hold for transport in a well-defined 
thermal environment. Without this condition one could 
still express the mobility by an integral over some correla- 
tion function as a minor modification of our proof would 
show. The physically transparent connection to the dis- 
persions emphasized here, however, would then be lost. 
Even though our relation is remarkably reminiscent to 
the well-known linear response result, a crucial difference 
should not go unnoticed. For a non-equilibrium steady 
state as investigated here, the relevant correlations in- 
volve a "sensitivity- weighted" current. As an observable, 
the latter requires knowledge of how the rates depend 
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